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Q,^ . Abstract. We report an ultra-high-resolution simulation that follows 

evolution from the earliest stages of galaxy formation through the pe- 
riod of dynamical relaxation. The bubble structures of gas revealed 
^ ■ in our simulation (< 3 x 10* years) resemble closely the high-redshift 

fT^ ' Lyman a emitters (LAEs). After 10^ years these bodies are dominated 

, by stellar continuum radiation and look like the Lyman break galax- 

04 ' ies (LBGs) known as the high-redshift star-forming galaxies at which 



point the abundance of elements appears to be solar. After 1.3 x 10^" 
years, these galaxies resemble present-day ellipticals. The comparisons 
of simulation results with the observations of elliptical galaxies allow 
us to conclude that LAEs and LBGs are infants of elliptical galaxies or 
bulge systems in the nearby universe. 



o 

:^ 

S 1 Introduction 

• Recent progress in observational devices and techniques has enhanced our knowl- 
J> [ edge of formation and evolution of galaxies on a firm statistical basis. Optical 
k>( " observations have revealed the presence of a number of LAEs at redshifts of z > 
rN ] 3 (Dey et al. 1998; Taniguchi et al. 2003 and references therein) as well as LBGs 
■ at redshifts of 3 < z < 6 (Steidel et al. 1996; Giavalisco 2002 and references 
' ' ' ' therein). A good fraction of LBGs exhibit Lya emission lines strong enough to 
be detected by narrow band observations (Shapley et al. 2003). Since LAEs and 
LBGs are quite young, they could hold direct interpretable information on the 
early chemical enrichment of galaxies, contrary to present-day galaxies which have 
undergone intense interstellar medium (ISM) recycling, thus erasing most of the 
early chemical history. 

The modern paradigm of the galaxy formation, based on the cold dark matter 
hypothesis, deduces that galaxies formed hierarchically in a bottom-up fashion. 
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where a larger system results from the assembly of smaller dark matter halos. 
Baryonic gas falls into the gravitational potential of dark matter halos, and con- 
denses rapidly as a result of the radiative cooling for atoms or molecules. The 
dense cooled gas clouds are the birth sites of stars, and massive stars born there 
explode as type II supernovae (SNe) in a few times 10'' years. Continual multiple 
SNe generate hot bubbles enriched with heavy elements and shocked dense shells. 
Hence, to study the chemical enrichment, it is crucial to resolve accurately the 
thermalization of the kinetic energy released by multiple SNe and the mixing of 
heavy elements. So far, the theoretical models of galactic chemical evolution have 
often assumed the homogeneous ISM (one- zone model), with the instantaneous 
and perfect mixing of heavy elements synthesized in SNe. However, the energy 
input and metal ejection by SNe are likely to proceed in an inhomogeneous fashion 
(Mori et al. 2004). Thus, simulations that can resolve SN remnants are required 
to properly model the chemical evolution of primordial galaxies. In this paper, 
we perform ultra-high resolution hydrodynamic simulations of a very large burst 
of multiple SN explosions in a forming galaxy. Also, by incorporating spcctro- 
photometric modeling with the simulations, the results can be directly compared 
to the observations on high-redshift forming galaxies. The outline of this paper 
is as follows. In Section 2, the numerical method is described. In Section 3, we 
present the simulation results and in Section 4 the possible link among LAEs, 
LBGs and elliptical galaxies is discussed. The full description of this study can be 
found in Mori & Umemura (2006). 

2 Numerical method 

Our simulation uses a hybrid iV-body/hydrodynamics code which is applicable to 
a complex system consisting of dark matter, stars and gas. The gas is allowed to 
form stars and is subject to physical processes such as the radiative cooling and 
the energy feedback from SNe. The gas is assumed to be optically thin and in 
coUisional-ionization equilibrium. Radiative cooling is included self-consistently 
with metallicity, using the metallicity-dependent cooling curves by Sutherland & 
Dopita (1993). The collisionless dynamics for dark matter particles and stars 
is treated by the A''-body method and the gas dynamics is pursued by a three- 
dimensional AUSM-DV scheme that can treat shocks with high accuracy (Mori 
et al. 2002). Since this scheme has a great advantage due to the reduction of 
numerical viscosity, fluid interfaces are sharply preserved and small-scale features 
can be resolved. 

Stars are assumed to form in rapidly cooling and Jeans unstable regions at a 
rate which is inversely proportional to the local dynamical time (see Mori et al. 
1997, 1999). When a star particle is formed, we identify this with approximately 
10^ single stars and distribute the associated mass of the star particle over the 
single stars according to Salpeter's initial mass function. The lower and upper 
mass limits are taken as 0.1 and 100 Mq, respectively. When a star particle is 
formed and identified with a stellar assemblage as described above, stars more 
massive than 8 Mq start to explode as Type II SNe with the explosion energy 
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of 10"''"'^ ergs and their outer layers are blown out with synthesized metals leaving 
the remnant of 1.4 Mq. Consequently, once a new star particle is formed, the 
energy, metals and material from Type II SNe are subsequently supplied to 8 cells 
surrounding SN region. We compute the chemical evolution using the calculations 
of stellar nucleosynthesis products by Tsujimoto et al. (1995). A mass of 2.4 Mq 
of oxygen is ejected from a Type II SN explosion. 

3 Simulation 

Following a standard cold dark matter cosmology with cosmological constant 
(ACDM), where we assume = 0.3, = 0.7, Vtb = 0.04, and a Hubble constant 
of Hq = 70 km s~^ Mpc~^, we consider the dynamical and chemical evolution of 
a protogalaxy with the total mass of 10^^ Mq. We assumed that the total mass of 
gaseous matter is 1.3 x 10^° Mq initially. If we suppose a 2a density fluctuation, 
it decouples from the cosmic expansion and begins to contract at redshift z = 7.8 
with the radius of 53.7 kpc. The angular momentum is provided by a uniform 
rotation characterized by a spin parameter of A = 0.05. Prior to this galaxy-scale 
fluctuation, subgalactic dark halos collapse and are virialized. According to the 
ACDM cosmology, twenty subgalactic condensations with mass of 5.0 x IO^Mq 
and radius of 8.6 kpc are distributed within the galaxy-scale fluctuation. The 
subgalactic virialized halos are assumed to follow the Navarro-Prenk- White den- 
sity profile (Navarro et al. 1997). The hydrodynamic processes are pursued with 
1024^ grid points. The simulation box has a physical size of 134 kpc and the 
spatial resolution is 0.131 kpc. 

Fig. 1 shows the result for the time sequence of the chemical enrichment, where 
the distributions of the logarithmic density, the temperature and the velocity, and 
oxygen abundance [0/H], are presented until 10^ years. In the first 10^ years, stars 
are formed in high-density peaks within subgalactic condensations and the burst 
of star formation starts. Then, massive stars in the star forming regions explode 
as SNe one after another, producing expanding hot bubbles surrounded by cooled 
dense shells. The gas in the vicinity of SNe is quickly enriched with ejected heavy 
elements, but a large amount of gas still retains low heavy-element abundance. 
Consequently, the metallicity distribution becomes highly inhomogeneous, where 
gas enriched as —5 < [O/H] < —1 coexists with virtually primordial gas. Since 
the density of the ISM is lower in the outer regions of subgalactic condensations, 
the expansion of hot bubbles is accelerated there and SN-driven shocks collide 
with each other to generate super-bubbles of ^ 50 kpc, and the surrounding high- 
density, cooled (T = 10"^ K) shells form at 3 x 10® years. The hot bubbles expand 
further by subsequent SN explosions, and the shells sweep up the partially en- 
riched ambient gas. The gas density in dense shells increases owing to the efficient 
radiative cooling mainly through coUisional excitation of neutral hydrogen. 

After 5 x 10^ years, the hot bubbles blow out into the intergalactic space. This 
SN-driven outflow is an efficient mechanism to enrich the intergalactic mediiim 
with heavy elements over a large cosmological volume (see Mori et al. 2002). On 
the other hand, the dense shells undergoes hydrodynamic instabilities induced by 
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Fig. 1. Spatial distributions of the stellar density, the gas density, and the oxygen abun- 
dance [0/H] of gas. The four panels in each column depict the time evolution of the 
simulation results until 1 Gyr. 



shell-shell interactions and radiative cooling, eventually fragmenting into cold fil- 
aments and blobs. These interactions are giving rise to an intricate multiphase 
structure in the inner halo, where 10^~^ K gas coexists with a cooler 10^ K phase 
from which it is separated by cooling interfaces. New stars are born in this enriched 
gas and again heavy elements are ejected from subsequent SNe. The rightmost 
panels show the structure at 10^ years. By this stage, the ISM is recycled repeat- 
edly and about 72% of the initial gas has been processed into stars. Eventually, 
some amounts of cool, dense filaments are left at the center. But, the most of vol- 
ume is filled with rarefied gas {n < 10~^ cm~^) that has intermediate temperature 
(lO**'^ K < T < 10^'^ K). At this epoch, the mixing of heavy elements is nearly 
completed. 



4 Discussion and conclusion 



The emission properties of the gas components are calculated for an optically- 
thin, coUisionally-ionized gas using MAPPINGIII code by Sutherland & Dopita 
(1993), and those of the stellar components are calculated using the evolutionary 
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Fig. 2. Predicted emission properties of tlie gas components assuming for an optically 
thin gas in collisional- ionization equilibrium (thin lines), and that of the stellar compo- 
nents using the evolutionary stellar population synthesis (thick lines). 



stellar population synthesis code PEGASE by Fioc & Rocca-Volmerange (1997). 
In practice, to obtain the SED, we sum up the SED of each grid point for the gas 
components and each star particle for the stellar components. Thus, the present 
simulations can be directly compared to the observations. 

Fig. 2 shows the spectral energy distribution (SED) of the simulated galaxy. 
Thin lines and thick lines depict, respectively, the gas emission and the stellar 
emission. The Lya emission comes mainly from high-density cooling shells, and its 
luminosity is more than 10^^ erg s~^ in the first 3 x 10^ years. This Lya luminosity 
perfectly matches that of observed LAEs (Matsuda et al. 2004). Furthermore, we 
find that the physical extent of ~ 100 kpc and the bubbly structure produced by 
multiple SNe are quite similar to the observed features in Lya surface brightness 
distribution of LAEs. After 3 x 10* years, the Lya luminosity quickly declines to 
several 10^^ erg s~^ that is lower than the observed level. This result suggests that 
LAEs can correspond to an early SN-dominated phase before 3 x 10* years. 

As seen in Fig. 2, after the sparkling phase of a primeval galaxy, the SED 
is dominated by stellar continuum emission, since the emission from cooling gas 
decreases immediately owing to the leak of explosion energy through the blowout 
of super-bubbles. The galaxy in this phase is featured with diffuse, asymmetric 
structures, and outflows of 100 ~ 500 km s~^. The total mass of long-lived stars 
is 9.3 X 1O^M0, and the mass of 1.5 x 1O^M0 is involved in the outflows at 2; = 3. 
These features look quite similar to those observed for LBGs (Adclbcrger et al. 
2003). In the light of such properties, the simulated post-starburst galaxy that 
has the age of 10^ years can correspond to LBGs. 

The; following dynamical evolution is studied by an TV-body simulation with one 
million particles, to explore the end-product of this galaxy. As a result, it is found 
that the assembly of sub-condensations and the virialization of the total system 
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arc almost completed in 3 x 10^ years, so that the system becomes in a quasi- 
equilibrium state. The resultant stellar system forms a spherical virialized system. 
The projected surface brightness distributions have a large central concentration 
that well accords with de Vaucouleurs' r^/"^ profile and the resultant absolute 
magnitude in blue band and visual band are Mb = —17.2 mag and My = —18.0 
mag, respectively. The colour U — V = 1.15 and V — K = 2.85 are consistent with 
the colour-magnitude relation of elliptical galaxies in Coma cluster of galaxies 
(Bower et al. 1992). Furthermore, the combination of the surface brightness, the 
effective radius = 3.97 kpc, and the central velocity dispersion ctq = 133 km s~^ 
is on the fundamental plane of elliptical galaxies within their scatters (Djorgovski 
& Davis 1987). These comparisons of the simulation results with the observations 
allow us to derive a significant conclusion that LAEs and LBGs are progenitors 
of present-day elliptical galaxies, and the on-going, major chemical enrichment 
phases. 
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